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Abstract 


This  external  client  report  summarizes  the  results  of  a  data  collection  campaign 
performed  with  the  CATSI  sensor  to  measure  the  detectability,  identification  and 
column  density  of  a  plume.  The  measurements  were  made  simultaneously  with  an 
imaging  sensor  known  as  FIRST,  which  was  developed  by  Telops  Inc.  The  objective 
of  this  measurement  campaign  was  to  provide  an  evaluation  of  the  Telops  FIRST 
spectrometer.  Measurements  were  made  at  the  500  m  open-air  corridor  at  DRDC 
Valcartier  on  September  29,  2004. 


Resume 


Ce  rapport  pour  un  client  exteme  rapporte  les  mesures  effectuees  avec  CATSI  sur  la 
detectabilite,  F identification  et  la  quantification  d’un  panache  de  SF6.  Les  mesures  ont 
ete  effectuees  simultanement  avec  le  spectrometre  imageur  FIRST  de  Telops.  Ces 
mesures  ont  ete  produites  afin  d’aider  revaluation  du  spectrometre  de  Telops.  Ces 
mesures  ont  ete  effectuees  sur  le  corridor  de  500  metres  de  RDDC  Valcartier  le  29 
septembre  2004. 
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Executive  summary 


The  development  of  techniques  for  the  passive  detection,  identification  and 
quantification  of  chemical  vapours  is  a  very  active  research  sector.  Telops  is  currently 
developing  an  imaging  Fourier- transform  infrared  spectrometer  for  this  application. 

This  external  client  report  presents  the  results  of  a  data  collection  campaign  held  on 
September  29,  2004  at  Defence  Research  and  Development  Canada  (DRDC)  - 
Valcartier.  The  results  consist  of  the  detection,  identification  and  the  integrated  path- 
concentrations  (ppm-m)  of  SF6  plumes  using  the  Compact  ATmospheric  Sounding 
Inteferometer  (CATSI)  and  the  Telops  imaging  spectrometer  (FIRST). 


Lavoie,  H.,  E.  Puckrin  and  J.-M. Theriault.  2005.  Passive  standoff  detection  of  SF6 
plumes  at  500  meters,  measurement  campaign  to  correlate  Telops  imaging  spectrometer 
(FIRST).  ECR  2004-373.  Defence  R&D  Canada  -  Valcartier. 
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Sommaire 


Le  developpement  de  techniques  pour  la  detection,  Fidentification  et  la  quantification 
en  retrait  de  vapeurs  chimiques  est  un  secteur  de  recherches  tres  actif.  Telops 
developpe  actuellement  un  spectrometre  imageur  infrarouge  par  transformee  de 
Fourrier  pour  cette  application. 

Le  present  rapport  pour  client  externe  presente  les  resultats  d’une  campagne  de  mesure 
effectuee  le  29  septembre  2004.  Les  mesures  de  detection,  identification  et  de 
quantification  fut  effectuees  sur  des  panaches  d’hexafluomre  de  souffre  situe  a  une 
distance  d’environ  500m.  Ces  mesures  on  ete  effectuees  simultanement  avec  les 
spectrometres  CATSI  de  RDDC  Valcartier  et  FIRST  de  la  compagnie  Telops. 


Lavoie,  H.,  E.  Puckrin  and  J.-M. Theriault.  2005.  Passive  standoff  detection  of  SF6 
plumes  at  500  meters,  measurement  campaign  to  correlate  Telops  imaging  spectrometer 
(FIRST).  ECR  2004-373.  R&D  pour  la  Defense  Canada-  Valcartier. 
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1.  Introduction 


This  external  client  report  (ECR)  presents  the  results  of  a  data  collection  campaign 
held  on  September  29,  2004  at  Defence  Research  and  Development  Canada  (DRDC)  - 
Valcartier.  The  results  consist  of  the  detection,  identification  and  the  integrated  path- 
concentrations  (ppm-m)  of  SF6  plumes  using  the  Compact  ATmospheric  Sounding 
Inteferometer  (CATSI).  These  results  will  serve  to  support  the  evaluation  of  the 
Telops  imaging  spectrometer  (FIRST). 

This  report  describes  the  monitoring  method  with  CATSI  and  the  parameters 
associated  with  the  chemical  plume  release.  This  is  followed  by  a  summary  of  the 
results  in  graphical  form  of  the  spectral  radiance,  column  density  (CL)  and  correlation 
coefficients  (r^). 
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2.  Material  and  methods 


2.1  Monitoring  method 

The  passive  standoff  monitoring  of  chemical  vapours  with  CATSI  is  based  on  the 
differential  detection  equation  derived  with  the  atmospheric  parameters  defined  in  Fig. 
1.  It  can  be  verified^  that  if  the  cloud  occupies  a  fraction / of  the  field-of-view  then  the 
radiance  differential  would  be  given  by, 

SL  j  =L  -Lj  =ALj  )(5  T  +N  -Lj  -ALj  ).  (1) 

Equation  1  represents  a  general  expression  that  represents  the  differential  radiance  of  a 
cloud  having  a  filling  factor  (/)  with  a  provision  for  handling  a  possible  background 
variation  {ALdear)  between  the  two  simultaneously  observed  scenes.  The  parameter 
Lgas  is  the  radiance  of  an  atmosphere  containing  the  cloud  where  Bgas  and  igas  represent 
the  Planck  radiance  and  transmittance,  respectively,  associated  with  the  cloud.  The 
expression  (1  -  denotes  the  cloud  spectral  emissivity,  and  Nmar  corresponds  to  the 
near- field  path  radiance.  For  the  general  case  of  a  gas  mixture  composed  of 
i  compounds,  igas  is  given  by. 


=  exp 


(2) 


where  ai(v)  is  the  spectrally  dependent  absorption  coefficient  (1/ppm-m)  of  compound 
/,  Ci  is  the  volume  concentration  (ppmv)  of  compound  /,  L  is  the  length  (m)  of  the 
cloud  sample,  and  v  is  the  wavenumber  (cm'^).  Note  that  the  quantity  of  interest,  C^L, 
is  given  here  in  ppm-m. 

The  basic  processing  strategy  of  the  detection  algorithm^  (GASEM  -  GASeous 
Emission  Monitoring)  consists  of  adjusting  a  set  of  gas  parameters  (e.g.,  CL  and 
temperature)  to  generate  a  calculated  spectrum  (SLcaic)  that  best  fits  the  spectral 
features  of  the  measured  spectrum  (dLmeas)  inside  a  spectral  band  from  v^in  to  Vmax-  The 
minimization  equation  can  be  written  as. 


’^max 


(3) 


Note  that  the  background  radiance  drift  between  the  two  scenes  (ALdear)  is  handled  by 
including  a  blackbody  temperature  variation  (A 9)  in  the  fitting  process.  In  the  current 
version  of  GASEM,  the  radiance  differential,  SLcaic,  is  calculated  using  Eq.  2  with  ai(v) 
taken  from  an  IR  spectral  database.  A  fast  atmospheric  transmission/emission  model 
has  been  integrated  into  GASEM.  This  model  is  a  look-up  table  version  of  the 
MODTRAN  modell 
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The  detection  and  identification  of  a  given  chemical  is  achieved  when  the  chemical 
cloud  parameters  retrieved  through  a  SIMPLEX  minimization  are  physically 
acceptable,  and  when  a  proper  correlation  factor  between  the  measured  spectrum  and 
the  spectrum  calculated  with  the  best-fit  parameters  is  obtained.  The  correlation  factor 
{/)  is  defined  here  as  the  square  of  the  Pearson’s  correlation  coefficient.  It  is  given 

by, 


=  - 


'^chan 

^(x,.  -x\yi  -y) 


i=\ 


^chan 


■Y.iyi-yf 


i=l 


i=l 


(4) 


ATMOSPHERE  BACKGROUND 


SENSOR  CLOUD 


Figure  1.  Diagram  and  the  associated  terminology  defining  the  geometry  for  the  differential  detection 

with  a  dual  beam  interferometer 


with  Xi  =  SLcaic,  yi  =  ^Lmeas,  ^nd  whcrc  richan  represents  the  number  of  spectral  channels 
of  the  chosen  band  for  the  spectral  fit,  X/  is  equal  to  the  spectral  radiance  differential 
calculated  with  the  best-fit  parameters  plus  the  background  temperature  drift,  AO,  and 
yi  is  equal  to  the  measured  spectral  radiance  differential.  The  terms  x  and  y  represent 
the  average  over  the  spectral  band  of  the  calculated  and  the  measured  radiance 
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differentials,  respectively.  This  correlation  factor,  which  varies  from  no  correlation  (r^ 
=  0)  to  perfect  correlation  (r^  =  1),  is  then  used  to  make  a  decision  on  the  presence  or 
the  absence  of  a  given  chemical  cloud.  Figure  1  shows  a  typical  picture  of  the  output 
screen  with  the  associated  parameters  generated  with  the  GASEM  software  during  an 
acquisition  at  a  spectral  resolution  of  8  cm  \ 

2.2  Chemical  plume  release 

In  order  to  obtain  a  relatively  uniform  vapour  cloud,  a  stabilized  source  of  SF6  from  a 
tank  was  connected  to  an  aluminium  pipe  with  a  diameter  of  4  cm  and  a  length  of  3  m. 
The  pipe  was  mounted  vertically  on  a  tripod  at  a  height  of  approximately  2  m.  The  top 
of  the  pipe  was  blocked  and  a  series  of  5 -mm  diameter  holes  were  linearly  distributed 
at  10  cm  intervals  along  the  pipe  to  act  as  point  sources.  This  arrangement  provided 
an  extended  linear  source  of  2  m  of  vapour  as  opposed  to  the  point  release  directly 
available  from  the  tank  outlet.  A  regulator  and  flow  meter  were  used  to  stabilize  the 
SF6  flow  rate  at  values  ranging  from  25  to  200  L/min.  A  stabilized  flow  rate  was 
maintained  for  1  min  for  each  of  the  measurements. 

2.3  Collected  temperature 

Table  1  presents  several  useful  temperatures  measured  during  the  data  collection 
campaign  of  September  29^^.  The  ambient  temperature  (Tamb)  and  background 
temperatures  (natural  (Tmat),  space  blanket  (T^^)  and  black  shingle  (T^^))  were 
recorded  for  reference  and  AT  determination. 


Table  1.  Temperature  of  air  (Tamb)  dnd  different  backgrounds:  natural  (Tnat),  reflective  (Tsb)  snd 

absorbent  (Tbs) 


TIME 

Tamb 

Tnat 

Tsb 

Tbs 

10h35 

14.1 

-36 

34 

15h05 

18.8 

17 

-36 

50 

16h05 

19.0 

17 

-36 

41 
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3.  Measurements of  SFe plumes at 500 meters 


The  spectral  signature  of  SF6  is  presented  in  Figure  2.  This  signature  consists  of  only 
one  band  at  947  cm'\ 

Figures  3  to  13  summarize  the  main  monitoring  results  for  the  spectral  radiance,  the 
column  density  (CL)  and  the  correlation  coefficients  (r^).  A  three-dimensional 
overlay  plot  for  each  single  spectrum  collected  during  the  SFs  release  was  generated 
for  each  measurement.  The  calculations  of  the  integrated  path-concentrations  and  the 
global  r^,  which  is  a  measure  of  the  quality  of  the  fit,  are  also  presented  for  each  gas 
release.  These  calculations  were  made  with  the  GASeous  Emission  Monitoring 
(GASEM)  algorithm.  Each  figure  is  also  tagged  with  the  acquisition  time  associated 
with  the  CATSI  sensor  (time  difference  between  FIRST  and  CATSI  measurements 
equal  150  s). 


Figure  2.  Spectral  signature  of  SFe 
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Figure  3.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATSI  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  100  L/min.  Measurement  performed  at  15h  16 

(CATSI  time)  on  natural  background. 
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Figure  4.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATSI  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  200  L/min.  Measurement  performed  at  15h27 

(CATSI  time)  on  natural  background. 
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Figure  5.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATSI  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  100  L/min.  Measurement  performed  at  15h34 

(CATSI  time)  on  natural  background. 
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Figure  6.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATSI  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  50  Umin.  Measurement  performed  at  15h38  (CATSI 

time)  on  natural  background. 
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Figure  7.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATSI  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  25  Umin.  Measurement  performed  at  15h41  (CATSI 

time)  on  natural  background. 
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Figure  8.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATSI  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  100  U  min.  Measurement  performed  at  15h56 

(CATSI  time)  on  reflective  background. 
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Figure  9.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATSI  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  50  Umin.  Measurement  performed  at  16h01  (CATSI 

time)  on  reflective  background. 
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Figure  10.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATSI  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  25  Umin.  Measurement  performed  at  16h05  (CATSI 

time)  on  reflective  background. 


DRDC  Valcartier  ECR  2004-373 


13 


X  10 


Figure  11.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATS  I  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  100  L/min.  Measurement  performed  at  16h  17 

(CATS I  time)  on  absorbent  background. 
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Figure  12.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATSI  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  50  Umin.  Measurement  performed  at  16h20  (CATSI 

time)  on  absorbent  background. 
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Figure  13.  3D  Overlay  plot  of  all  radiance  spectra  (Watts/cm^-sr-cm'^ )  collected  with  CATSI  sensor  (top 
view),  graph  of  the  related  CL  (blue  curve)  and  global  (red  curve)  (lower  graph)  for  a  measurement 
performed  on  a  plume  generated  with  a  flow  rate  of  25  Umin.  Measurement  performed  at  16h26  (CATSI 

time)  on  absorbent  background. 
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4.  Conclusion 


This  external  client  report  presented  the  spectral  measurements  performed  with  CATSI 
on  September  29,  2004  for  the  release  of  SF6  at  500  m.  The  measurements  were 
obtained  simultaneously  with  the  Telops  imaging  sensor  (FIRST)  to  validate  the 
capability  of  the  sensor  and  to  acquire  a  data  set  for  the  development  and  validation  of 
the  detection,  identification  and  quantification  algorithms. 

The  Gasem  algorithm  was  used  to  detect,  identify  and  quantify  the  simulant  (SF6) 
plume  released.  The  algorithm  was  setup  to  fit  the  column  density  (CL)  and  the 
spectral  slope  in  function  of  the  wavenumber  (ATa).  The  column  density  was  obtained 
for  different  background  type  and  different  simulant  flow  rate. 

The  column  densities  (CL)  obtained  for  a  released  plume  over  natural  backgrounds 
were  consistent  with  values  obtained  in  other  field  trials.  The  column  densities 
obtained  for  the  artificial  backgrounds  were  inconsistent  with  the  values  obtained  for 
natural  backgrounds,  being  as  high  as  an  order  of  magnitude  greater  than  expected.  At 
the  moment,  there  is  no  obvious  explanation  for  this  difference. 
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